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Abstract 

Validity of the coalescence model of exotic atom formation in the ion col- 
lisions is studied . An estimate is given for the rates of P state production in 
the electric fields of the colliding ions . Possibilities of quantum beat effects are 
indicated . 

1 Introduction 

We study possibilities of production of pairs composed of an elementary particle 
and its antiparticle in the atomic bound states. The particles may be produced in 
proton or ion collisions and the simplest cases of practical interest are the lightest 
pairs e+ e~ (positonium) , 7r+ tt~ (pionium) or K~ (kaonium). These atoms 
are large systems when compared to the regions where the pairs are created . 
Thus , amplitudes for the formation of atoms in final states are proportional 
to the atomic wave functions at the origin and the pair creation amplitudes . 
This approximation known as the coalescence model yields rates of mesic atom 
formation high enough to allow experiments . Thus the search for pionium has 
been undertaken or contemplated at CELSIUS, CERN, COSY, Indiana UCL, 
[1,2,3,4,5]. 

The main physical motivation to study pionium (or kaonium) is a mea- 
surement of its lifetime . In the S states the latter may be expressed by the 
appropriate particle-antiparticle scattering length [6] , a quantity of great the- 
oretical interest . The main difficulty of these experiments is that the decay 
of the IS state , due to strong interactions , is very fast ( about 10"^'^ sec for 
pionium) . In order to find the lifetime one has to compare the main tt'^ tt^ decay 
mode with very rare two gamma decay mode . In principle ,the atom's life may 
be prolonged by orders of magnitude if it is formed in the P states that decay 
slowly to the ground state by the X-ray emission [7] . This would open new ex- 
perimental possibilities and a new scale given by the calculable X-ray decay rate 
. However, a direct P wave pair production with the subsequent P state atom 
formation in a coalescent manner is very unlikely [6] . In this paper another 
mechanism for the P wave atom production is considered . We assume the pair 
to be produced at short ranges in an S wave . Subsequent formation of the P 
state atoms is achieved by a Stark mixing in the external electric field of the 
colliding particles . A related problem discussed here is the rate of destruction 
of the S state atoms formed via the coalescence mechanism and destroyed by 
the external electric field . The second question is relevant to the studies of 
atomic production in the proton-nucleus collisions at high energies [4] . 
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The atom formation is discussed here as a problem of a pair created in an 
external static Coulomb field . In the full complication it becomes a formidable 
problem of three body Coulomb interactions in the continuum. Only part of it 
has been solved with satisfactory precision. In particular we consider only atoms 
moving with high velocities. In a kind of quasi-clasical solution we find the rate 
of destruction of the IS atoms in the external field and the rate of production of 
the P wave states. The former is significant the latter is rather small. An upper 
limit for these quantities is given with a help of completeness approximation . 
A system of succesive approximations is shown to converge to this limit. 

This model shows that the participation of different atomic states depends 
strongly on the final momentum of created pairs. If only a limited number 
of final states is allowed the atomic system displays effects of quantum beat . 
The atomic pair in the external, polarising, electric field forms quasi-stationary 
Stark type combinations. As a consequence the probability of production of 
each atomic state oscillates with respect to the final atom momentum. This 
phenomenon is rather unstable against the approximations. We believe ,it de- 
serves more extended study as it may be related to the production of correlated 
e+ e~ pairs seen at GSI ,[8] , in the heavy ion collisions. 



2 The model of atomic formation 

This section presents basic assumptions of this model. The equations that de- 
scribe the atom formation are derived and simpified solutions are found. The 
latter allow to form an intuitive picture of the formation mechanism. An ap- 
proximation scheme for the high energy region is derived. 

In cases of practical interest the atomic Bohr radii , B , are much larger 
than radii of the interaction region where the pair is created. This leads to 
the expectation that the atom formation amplitude A(E) is proportional to the 
wave function at the origin : 



A{E) = F^p{e = 0,f=0) (1) 

where F is an amplitude for the pair creation . 

This paper is devoted entirely to the determination of the final state wave 
function tp. The latter depends on two coordinates : ^that describes the center 
of mass of the pair tt"*" tt" , and r that is the relative coordinate in the center 
of mass. We assume the adiabatic situation i.e. fixed ions and fast mesons . 
Then the external polarization potential V of colliding nuclei acting on the pair 
becomes : 

V = Za{—^-—^) (2) 
where Z is the total charge of the nuclei. 
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Let us start with the Lippman-Schwinger integral equation describing non- 
relativistic dynamical formation of atomic states of the pair tt+tt". We have, 



^ = ^0 + j GV^d-'gd-'r (3) 

and the boundary conditions , corresponding to the ingoing waves , arc to be 
specified later .In this equation the "internal " interactions within the pair have 
been solved and described by the Green's function G . At first only the discrete 
terms of this Green's function are taken into account. 

G = G{r, g, ?, g')f^^\ (pnim >< Vnim I G{g, g', E - En), (4) 

nlm 

where E is the total energy of the pair . We are interested in the high kinetic 
energies ( MeV region ) and neglect the atomic binding energies En ( KeV region 

2 

). Then , E = j^. where m is the paticle mass. In addition, we assume that the 
pair is created with zero total angular momentum. Hence, the wave function 
can be written in the following form: 

nlm 

and the radial wave functions are assumed in the form 

^nl 'Pnl{r)(f)nl{g) (6) 

(t>ni = jiipg)ani{g) (7) 

where the unknown functions ani{g) are to be found . The form of this expansion 
although general has been chosen for the high-energy region which allows simple 
approximate solutions for the functions a . After eliminating the internal degrees 
of freedom by integration over r one obtains an infinite set of second order linear 
differential equations for a„i that involve effective potentials in the "external" 
variable g 



where 
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Delta{l,l',X) = {2l + l){2l' +l){2X+l)i^[ J (9) 

and f p Q Q j is the Wigner symbol. 

These potentials are related to the external Coulomb potential V by the 
multipole expansion : 
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V{r, Q)=Y^ VyX;^{r)Yxm{Q) (10) 

Am 

with V\ the radial part of the A-pole polarisation potential given by 



l^.(r,.) = ^^,A = l,3... (11) 



where x is g or | . 

The dominant dipole effective potentials are long-ranged of the q~'^ asymp- 
totics . At distances close to the atomic radii these potentials fall down sharply 
and exhibit zeros of high order at the origin . Few examples are given in Ap- 
pendix A. The long range of these potentials allows some vital approximations 
in the high energy region . We refer to Appendix B where the details are dis- 
cussed in the language of integral equations while some numerical studies are 
performed in the next section .One effect of the range is some smoothness of the 
functions a and this allows for differential equations of the first order. Another 
special feature is related both to the long-range and the angular momentum car- 
ried by the effective potential . This allows for the expansion (6,7) where , one 
would otherwise expect also some hankel functions corresponding to ingoing( 
outgoing ) waves . The first order equation discussed in Appendix B is 

<i = 51 |-^";"'''«"'H-iy' (12) 

where the length is expressed in the Bohr radius units . These equations for 
o-nlio) can be written in the following matrix form 

a' = Ua (13) 

or 

ia' = iUa (14) 

where U is an antisymmetric matrix and respectively iU is hermitian. The 
boundary conditions can be formulated via integral equation 



a = Oo -h 



/"OO 

/ Ua dg (15) 

•I o 



which is the set of integral equations derived in Appendix B . For example for 
the incoming wave coresponding with the state IS 

ao = (1,0,0,...) (16) 

and the vector a represents following descrete states 

a = [IS, 2S, 2P, 35, 3P, 3D, ...) (17) 
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3 Illustrative example - two discrete states 



To illustrate the formalism introduced above let us consider an idealised atomic 
system of only two: 2S and 2P states . It is easy to find ,m this case , the exact 
as well as the approximate solutions for the integral or differential equations for 
the final atomic state of the pair . If the incident wave is formed with the 2S 
state one obtains from eq(12) or eq(15) : 



a2o{g) = cos{— / V2120) (18) 
2p Je 

which generates the S wave amplitude of the pair with respect to the formation 
region . Now the 2S atom formation probablity pr(2S) is given by the square of 
the amplitude at the origin 

Z f°° 

pr{2S)=Acos''{—j ^2120) (19) 

where A is a constant involving pair production amplitude and some phase space 
integrations , [ see 6 for more details]. 

In the same way one can solve for an amplitude a2o in the case of the incident 
2P state 

Z f°° 

a2o{g) = sin{— J ¥2120) (20) 

and the probability of the 2P atomic state formation is 

Z f°° 

pr(2P) = Asin^— / ^2120) (21) 

The total probability of the atom formation is then energy independent 

pr = pr{2S) + pr{2P) = A (22) 

and equals in the high energy limit to the 2S state formation probability given 
by the coalescence model . In figs 1,2 parts of the plots of pr(2S) and pr(2P) are 
presented for a special case of pozitonium formed in the strong field ( Z=170) 
of heavy ions . 
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nonrelativistic discrete 2p-2s model, creation of 2s 
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Let us notice several essential points. 

1) The oscillations of the atomic production probability reflect quantum 
beating of this system . In the strong external electric field there are two di- 
agonal combinations S(i)S(e)+P(i)P(e) and S(i)S(e)-P(i)P(e) where (i) denotes 
the internal atomic state and (e) the state of the external angular momentum 
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of the atom . These two components beat in time and in space when the atom 
propagates through the polarising electric field. 

2) This simple case may be solved exactly . Wc find that our high energy 
approximation works well for kinetic energies higher than the barrier inherent 
in the effective potential . In practice it covers the region of three last peaks in 
the formation probability . 

3) The beating effects persist at kinetic energies higher than the atomic 
binding energies by orders of magnitude . The Coulomb potential and the 
dominant dipole polarising potential have no inherent scale parameter. This 
parameter is introduced via the atomic matrix elements (8) that determine the 
effective V . Thus the energy scale in fig 1 is determined by the dimensionless 
parameter Z/pB . It is much larger than the scale of the atomic binding energies 

The question to be answered is to what extent this beating of the Stark mixed 
atomic state persists in complete calculations . Such calculations should allow 
decays of the incident state to higher states and to the continuum . We study 
it by expanding the basis of the discrete levels allowed and by an approximate 
treatment of the continuum .The latter is discussed in the next section . 

4 Approximate treatment of continuum 

In appendix B we estimate the continuum part ip'^ of the wave function ?/' • 
This part is related to the continuous spectrum involved in the Coulomb Green 
function entering the Lippman-Shwinger eq(3) . We argue that at high 
energies there exist an approximate relation 

< ^nlm I VCV I ^Pn'l'm' > « < ^Pnlm \ V5{f- ?)GpV \ Ifin'l'm' > (23) 

Then we expand the delta function into a system of eigenstates 

d{f- r') = ^ I i>nim. >< ipnim | +conUnuous part (24) 

nlm 

and neglect the continuous terms to obtain the following approximation 

< Vnlm I VCV I ipn'Vm' > « < V>nlm | VC^V \ <p„>l'm' > (25) 

where is the discrete part of Coulomb Green's function. In this way , we 

expand the continuous piece of the wave function ip'^ into discrete states in the 
same way as it was done in eqs (5,6,7) with the ip"^ . The coefficients of this 
expansion are named b . Now defining a new vector o 

a = (a, b) (26) 

we obtain the equation 

a' = Ua (27) 
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where the new potential U is defined as follows: 

The equation for a is soluble with simple Runge-Kutta methods. This procedure 
doubles the dimension of the previously described set of differential equations 
.Its effect on the discrete part is equivalent to an optical potential contribu- 
tion leading to the equation 

ip"^ = ijf^ + C^Vi/ + G'^VG'^ViIj'^ (29) 

It turns out that contributions from ip'^ at the origin vahish quickly with in- 
creasing energy and for intermediate energies stop plaing an important role. 



5 Calculations and Conclusions 

This calculation decribes formation of atomic states of a pair e+ e~ , 7r+ tt" etc. 

in the electric field of colliding nuclei that create the pair. It is shown that the 
rate of creation of pairs in a given atomic state depends strongly on the final 
energy of the pair. This result depends on one scale parameter so it can 

be rescalcd for different masses of the pair, total atomic numbers of colliding 
nuclei Z and energies. The probability to form the IS state changes rather 
smoothly with energy from negligible values at low energies to the coalescence 
model result at high energies . A reliable calculation of the formation rates in 
higher atomic states requires a large number of atomic states in the basis of 
expansion . We have solved the system of differential equations up to K= 21 
states. Good convergence of probabilities with respect to K have been obtained 
only for the first three states. Within higher shells some traces of the beating 
effects are still visible although these are rather unstable with respect to K . 
This question requires further studies and possibly better approximations as for 
large K the effective potentials create numerical difficulties . 
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6 Appendix A - Wave functions and potentials 



The first three wave functions and effective potentials are given below . The 
atomic units ( B =1) are used. The radial wave functions : 

^10 = 2e-'- (30) 

(31) 

v., = (32) 

generate , under the dipole external field , the following transition potentials : 



, , (5.959 +(-5.959- 17.88a; - 26.82a;2-20.11a;3)e-3=^ 
^2110 (a;) = ^2 (^^) 

, , (-24. + (24. + 48.a; + 48.a;2 + 32.a;3 + 16.a;^)e-2^ 
V"2i2o = ^ 2 34 

, , (4.33 + (-4.33- 7.22a; -6.02a;2-3.17a;3-l.lla;^-3.71a;5)e-5^ 
^3021 (a;) = 

(35) 

All other transition potentials defined in section 1 can be found in explicit 
form as elementary functions. 



7 Appendix B - High energy approximations 

The wave function V' can be divided into two parts - discrete and continuous: 

VJ = v/ + V (36) 

This separation follows the same seperation of the Coulomb Green function. 

G = G'^ + G'' (37) 

Where the first part contains the atomic and the second part the continuous 
internal states of the pair . Now Lippman-Schwinger equation is equivalent to 
coupled equations for the two parts oioitp. 

In this calculation we are interested in atomic states , the continuum will be 
described in a simplified way similar to the optical potential method . Thus we 
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neglect the G'^Vtj^'^ term that changes the continuum part of the propagator . 
Now 

■4,'' = C'VtP^ (39) 
is substituted into equation for tjj" that generates the optical model equation 
iP'^ = 'tP^ + G'^Vi)'^ + G'^VCVi}'^ (40) 

For high energies, much larger than the atomic binding energies , one has an 
approximate relation for the matrix element 

> (41) 

so in the above equations wc interchange G'^ with Gq. The discrete part of 
propagator has the following form 

G'^ = G'^ir, g, r', g') = J2\ "P^im >< ^nim \ Gp{g, g', E - En), (42) 

nlm 

where for high energies E we neglect the binding energies E^- In a similar 
fashion for small relative momenta of the pair the Green function Gq will be 
used in the following approximate form 

Go(r, g, P, g') = 5{r- r')Gpig, g'), (43) 

The external potential changes the internal angular momentum . This is 
crucial for our problem and allows an approximation that we discuss now . Let us 
go over to the angular momentum representation and follow the expansion given 
in eqs(5) and (6) . The propagator is also expanded into spherical harmonics 

GpiQ, g') ^JZaiis^ e')Yil{mm{g') (44) 

Im 

with 

9i{q, q') = "ip.n{PQ<)hi{g>) (45) 

Now we can obtain the system of integral equations for the radial functions </> 
of eg(6) 

'Pnlio) = jlo{PQ)5nno5llo + Z i gi{g, g') Vnln'l'4'nl 

■' n'l' 

+Z^l gi{g, e') J2 W^n"n'r(e', Q")9c.{g', g")<l^n'i' {g") (46) 

n'l'a 

The transition potential Vninn' is defined in the following way: 

Vnln'A0)= y < I I > nMn(l l' X} (47) 
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and W^i^n, is given by the formula 

^ < ^ni I Vx{e,r)VA9' .r) I ^u'v > Delta{l,a,X)Delta{l\a,X) (48) 
A,Atr,3... i2a + lW{2l + l){2l' + l) 

Iterations of these integral equations involve products of the bessel and neu- 
man functions . These are integrated with the effective potentials V which are 
predominantly long-ranged . For high energies i.e. for high momenta p these 
functions oscillate rapidly and one obtains approximate relations for mean val- 
ues: 

<Mx)jv{x)>^^cos{{l-l')^) (49) 

1 TT 

< ni{x)ni>{x) ^cos{{l - I')-) (50) 

< ji{x)ni'{x) >« ^sin{{l - r)|) (51) 

where 

X = pQ (52) 

and symbols < f{x) > denote the mean value of an oscillating function f(x) 
weihted by V. We repeat again , that the success of these approximation is 
related to the long range of the dipole polarising potentials V . When we write 
(j)ni in the standard form 

^ni = 31 {PQ)ani {q) + hi {pg)bni {g) (53) 
and apply the expressions for mean values of bessel functions we obtain 

+ / -^siniia - l)^)sin{il' - a)|) g")n'ng")dg'dg" (54) 

^ n'l'ct 

The equation for 6„; has the similar form: 

Ki{g) = / :^sm((«-r)-) VKin';'n/(£'')<^e' 

+ /' - l)\>in{{l' - a)\) J2 W^in'i'iQ', Q")bn'v{g")dg' dg" (55) 

but because there is no initial term 5nno^m the solution of the above equation 
is trivial - i.e. 6„; = 0. The equation for a„; is the gateway for the high energy 
approximation. 
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